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Abstract 

One interpretation of proton stability is that it implies the existence of extra-flat directions of the 
minimal supersymmetric standard model, in particular u c u c d c e c and QQQL, where the operators 
lifting the potential are suppressed by a mass scale A which is much larger than the Planck mass, 
A ^ 10 26 GeV. Using D-term hybrid inflation as an example, we show that such flat directions 
can serve as the inflaton in supersymmetric inflation models. The resulting model is a minimal 
version of D-term inflation which requires the smallest number of additional fields. In the case 
where Q-balls form from the extra-flat direction condensate after inflation, successful Affleck-Dine 
baryogenesis is possible if the suppression mass scale is ^ 10 31 — 10 35 GeV. In this case the 
reheating temperature from Q-ball decay is in the range 3 — 100 GeV, while observable baryon 
isocurvature perturbations and non-thermal dark matter are possible. In the case of extra-flat 
directions with a large t squark component, there is no Q-ball formation and reheating is via 
conventional condensate decay. In this case the reheating temperature is in the range 1 — 100 TeV, 
naturally evading thermal gravitino overproduction while allowing sphaleron erasure of any large 
B — L asymmetry. 
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I. INTRODUCTION 



Successful models of supersymmetric (SUSY) inflation should ideally satisfy a number 
of requirements: natural compatibility with supergravity (SUGRA), lack of fine-tuned cou- 
plings, successful post-inflation era including reheating and baryogenesis, and compatibility 
with unified models of particle physics. With respect to these conditions, SUSY hybrid 
inflation models have a particular attraction [J They can achieve sufficient inflation 
without requiring very small or fine-tuned couplings, and in the case of D-term hybrid in- 
flation they are naturally compatible with SUGRA [3]. Focusing on the D-term inflation 
case, a natural question is the origin of the fields in the D-term inflation sector. The U(l) 
gauge field and charged vector pair Q± of D-term inflation might be understood as compo- 
nents of an extended gauge theory. However, the inflaton is usually a gauge singlet which 
is added to the model for no other reason 1 . If we do not add such a singlet, can D-term 
hybrid inflation still occur? Here we argue that it can. The vector pair will naturally couple 
to any gauge-invariant combination of fields in the MSSM. Such gauge-invariant products 
(monomials) also characterise flat directions of the MSSM. Thus a natural possibility is that 
a flat direction can play the role of the inflaton in D-term inflation models 2 . In this model 
the number of additional fields required for inflation is reduced to just a U(l) gauge field 
and the Q± vector pair, so providing a minimal version of D-term inflation. As we will 
show, conventional MSSM flat directions lifted by Planck scale-suppressed gauge-invariant 
superpotential terms are unsuitable. This is because such terms generally lift the flat di- 
rection scalar at field strengths well below the value required for inflation. However, it is 
known that certain gauge-invariant superpotential terms must be suppressed by more than 
the Planck scale or forbidden entirely. The d = 4 operators u c u c d c e c and QQQL will lead to 
rapid proton decay if they are only Planck scale-suppressed [9(. One way this problem can be 
solved is by assuming that the underlying complete theory introduces a dynamical suppres- 
sion factor into the non-renormalisable superpotential interactions, such that the effective 
mass scale suppressing the dangerous operators is A > 10 26 GeV [9[. It is also possible that 



1 Models exist which attempt to identify the inflaton with a known field, such as a right-handed sneutrino 

flflB. 

2 An interesting model using MSSM flat directions^ inflatons, which has a quite different philosophy with 
respect to fine-tunings, is given in 0]. See also 
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this dynamical suppression will also apply to all higher-order MSSM superpotential terms 
lifting the flat direction, such as (u c u c d c e c ) 2 and (QQQL) 2 . We will refer to a flat direction 
for which this is true as an 'extra-flat direction'. An alternative interpretation of the absence 
of proton decay is in terms of a discrete symmetry which eliminates the dangerous d — 4 op- 
erators lo| . In this case it is possible that the higher-order operators will be unsuppressed. 
However, as we will show, such unsuppressed flat directions, even if higher-order, cannot 
serve as an inflaton. If the existence of extra-flat directions is the correct interpretation of 
the absence of proton decay in the MSSM, then an extra-flat direction scalar could serve 
as the inflaton in a D-term inflation model. The extra-flat direction potential at large field 
values is naturally lifted to an inflationary plateau by its gauge-invariant superpotential 
coupling to $ + <l>_. Reheating and possibly baryogenesis would then come from the decay of 
the flat direction inflaton, via either Q-ball decay or conventional homogeneous condensate 
decay, depending on the t squark component of the flat direction. In this paper we will 
study D-term inflation along an extra-flat direction of the MSSM. The paper is organised 
as follows. In Section 2 we discuss extra-flat directions and the resulting D-term inflation 
model. In Section 3 we discuss reheating and baryogenesis. In Section 4 we present our 
conclusions. 

II. D-TERM INFLATON ALONG EXTRA-FLAT DIRECTIONS 
A. Potential 

We consider a flat direction $ in the MSSM and introduce two additional fields $± 
charged under a U(l) gauge group with the Fayet-Illiopoulos term £. The superpotential is 

w = + 2 m 
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where Ai j2 are Yukawa couplings and M is the reduced Planck mass, M = M P i/\/8tt 3 . We 
will present results for general m and n, specialising to the case of most interest m = n = 4, 
corresponding to $ 4 ~ u c u c d c e c or QQQL. Proton stability in the case of u c u c d c e c or 
QQQL requires that Ai ^ 10 -8 , corresponding to an effective suppression mass scale A = 
Mj Ai ^ 10 26 GeV. However, A2 is unconstrained by phenomenology and will be determined 
by the inflation model. The scalar potential in the global SUSY limit is then 



V = 



A 2 n 



nM n ~ x 
9 2 



(I0 + | 2 +I0-| 2 ) + 
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The supersymmetric global minimum is located at 4 

(0,0 + ,|0_|) = (O,O, v / e). 



(2) 



(3) 



If 



AiA 2 m+ "- 2 



« m 2 



(4) 



is satisfied, the mixing between + and 0_ is negligible. The potential is then simplified to 
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A 2 r ' 



The critical value of is given by 
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fy(e+l0 + l 2 -^l 2 ) 2 . 



1/n 



(5) 



|A 2 



(6) 



which determines the stability of the 0_ field at the origin. The origin is a false vacuum for 
|0| > |0 C |, while it is unstable for \<f>\ < \<f> c \. 



3 A SUSY mass term W D ^.$+$_ has not been included. This term would induce n minima with 
nonvanishing VEV for which consist of squark and/or slepton VEV and lead to large baryon or lepton 
number violation in the MSSM. Although in most cases there is no symmetry which can exclude such a 
term, we note that for the case m = n this term can be excluded by an R-symmetry which allows the 
terms $™ and $™$ + <f>_. In addition, if /j, is less than the scale of soft SUSY breaking terms, the minimum 
of the potential can be at $ = 0, while for larger fj, there can be directions in the complex <!> plane along 
which the field evolution can avoid the minima with $ 7^ 0. 

4 Note that there is a SUSY flat direction when \i — and <j> = 0, such that |^+| 2 — \4>-\ 2 = £• However, 
the minimum with (f) + = is selected since cf) + gains a large mass when <fi ^ during inflation. 
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B. Inflationary expansion 



For |0| > 1 C |, 0_ = is a local minimum and there is the false vacuum energy from the 
.D-term, which drives inflation. The potential during inflation is given as 

V~-q 2 e ( 1 + ^-ln — ^ , (7) 

where a = \f2Re{(f)) is the canonically normalised inflaton and A* is the renormalisation 
scale. Inflation ends when the inflaton reaches the larger of a c = \/2\(j) rmc \ and 

_ VngM 

where <T/ corresponds to the end of slow-roll. However, a non- vanishing F-term potential is 
also present in this model. Hence, we need to ensure that the condition Vp <C Vd is satisfied, 
which requires that 



A,o"' 1 



« \g 2 e (9) 



Mm-3 

is satisfied. Note that when this is satisfied, equation is also satisfied. The dynamics of 
the inflaton field is similar to that in the minimal D-term hybrid inflation model The 
solution of the slow-roll field equations is 

W = * + do) 

Here, a Q = max[cr/,cr c ] is the expectation value of inflaton when inflation terminates. The 
spectral index is 

1 / 27T 2 cr n 2 V 1 

while the value of £ x / 2 normalised to the curvature perturbation is 

° / i + r-^7 • (12) 



M \ N J V ng 2 M 2 N J 
In the case of ctq = cr c , which we will show is true in examples of interest, we find 



27rv 2 __(\ 2 y /n 



ng 2 M 2 N 

with 



(13) 



^) (^1'" ( 14 ) 
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^2 = ^2 C corresponds to <7n = 2a c , with a ~ cr c throughout inflation when A2 < \2 C - When 
A2 ^> A2 C as well as the case of cxo = cr c , the spectral index is n s = 1 — 1/N ~ 0.98, as in 
conventional .D-term inflation, while the value of £ 1//2 required to account for the observed 
curvature perturbation (P^ 2 = 4.8 x 10" 5 ) is £ 1//2 = 7.9 x lO 15 ™ 1 / 4 GeV. On the other hand, 
in the case where A2 <C A2 C , the spectral index approaches n s = 1 while the value of £ 1//2 is 
reduced by a factor {X-z/ \2 c ) 1 ^ 2n ■ 



C. Comparison with observations 



The spectral index observed by WMAP, n s = 0.958 ± 0.016 (1-cr) [ll|, is substantially 
smaller the D-term inflation value. In addition, WMAP data permits at most an O(10)% 

nnn 

contribution to the CMB power spectrum from cosmic strings |12l . 1131 . [14J , which implies 



that £ 1/2 < 4 x 10 15 GeV. (Here we have used Gu = 2 x 10~ 6 for the I = 10 WMAP 

n 

normalised string tension |15l|.) One way to interpret the WMAP observations is that they 
correspond to an adiabatic curvature perturbation with n s w 1 combined with a 10% cosmic 
string contribution [16j, which can be achieved by making A2 sufficiently small compared 
with A2 C . In this case the apparent spectral index of the combined perturbation is effectively 
lowered and can be in agreement with the 3-year WMAP data analysis jl6|. It is a striking 
feature of D-term inflation models in general that they have a solution which increases n s 
while decreasing the cosmic string contribution, just as required for this interpretation of the 
WMAP observations. With respect to this possibility, the extra-flat direction model has a 
possible advantage over conventional D-term inflation. The contribution of cosmic strings to 
the CMB power spectrum is proportional to /i 2 = (27r£) 2 . In the case of conventional D-term 
inflation with n s — 1, the value of £ 2 in the limit A 2 A2 C is proportional to A 2 .. Therefore 
A2 must lie within a rather narrow range of values for the cosmic string contribution to be 
O(10)%. In the case of the extra-flat direction inflaton, the dependence is oc A^ n . Therefore 
the CMB contribution varies much more gradually with A2 e.g. £ 2 oc A^ 2 for the case n = 4. 
Thus an O(10)% contribution is obtained for a much wider range of A2, making it perhaps 
a more natural possibility than in conventional D-term inflation. 
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D. Constraints from cosmic string bound, SUGRA and potential flatness 



We first check that the 10% cosmic string condition £ x / 2 »s 4x 10 15 GeV can be satisfied for 
reasonable values of g when \<ft c \ 2 is small enough compared with M 2 for SUGRA corrections 
to be neglected. We will require that \<p c \ < kM, with k < 0.3, so that \<p c \ 2 < 0.1M 2 . From 
equation ([6]), this implies that 

< A 2 A:"M 

9 ~l&r- (15) 

For the case n = 4 and £ x / 2 k4x 10 15 GeV, equation ffl~5j) implies that 



^•^-j . (16) 

Thus A2 should not be small compared with 1 if g is not very small compared with 1. From 

equation (TT2|) . to suppress £ x / 2 from 7.9 x 10 15 GeV to 4 x 10 15 GeV in the case n = 4 we 
require that A 2c /A 2 ~ 250, which implies that 

A 2 « 8 x 10~V 3 . (17) 

Equations ffl~6]) and f|T7|) imply that 

p< 0.03 f-^ (18) 



0.3 



0.3 X 3 



A 2 > 0.03 ( — ] . (19) 



and 



Thus k > 0.1 is necessary when A 2 < 1 in order to satisfy equation ( fi~9l) . 0.1 < ^ 0.3 
then implies that g ~ 0.01 — 0.03 and 0.03 < A 2 < 1. Therefore, as in conventional D-term 
inflation in the small coupling limit, g must be somewhat smaller than the Standard Model 
gauge couplings [kJ. In addition, A 2 must be much larger than Ai 5 . We next evaluate 
Eq. to find the condition on Ai for F-term corrections not to spoil the flatness of the 
inflaton potential. In general we find 

A ' « • (20 » 



We have assumed that a c > at. For the case n = 4 this requires that e 1/2 /M > |A 2 | 5 3 /^ 4 - With 
£ x /2 ~ 4 x 10 15 GeV and g w 0.02 this is easily satisfied. 
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For the case m = 4 this gives 

Al « 7xl0 - 9 (^) 3 ( J - 1 ^)\ (21) 

Thus for values of Ai which satisfy the proton decay constraint, Ai < 10~ 8 , the flat direction 
potential is easily sufficiently flat to serve as an inflaton. However, for the case of unsup- 
pressed n = m = 4 flat directions with Ai ~ 1, the F-term would violate the flatness of the 
flat-direction inflaton potential. An alternative solution of the proton decay problem is to 
consider elimination of the m = n = 4 operators entirely by a symmetry. In this case we 
expect to have unsuppressed operators with m = n = 8, such that Ai ~ A 2 ~ 1. However, 
in this case the F-flatness condition will still be violated. For n — 8, equation (TT51) implies 
that 

/ k \ 8 

g < 5 x 10- 3 A 2 f — J . (22) 
To suppress £}l 2 to 4 x 10 15 GeV, with n = 8 we need A 2c /A 2 ~ 6.5 x 10 4 . This implies that 

A 2 = 9 x lO^V 7 • (23) 
Therefore if A 2 < 1 we have g ^ 0.03. Equation fl23|) combined with equation (1221) gives 

4~°- 02 fAl • ( 24 ) 



s 0.3 / 

Thus k ^ 0.5 is necessary if g ^ 0.03. The F-term flatness condition equation ([20]) for m = 8 
is 

^-- 6 («4)(t)'(^v) • (25) 

Thus the m = n = 8 flat direction will also need to be extra-suppressed to have a flat 
inflaton potential, even if the m = n = 4 term is completely eliminated by a discrete 
symmetry. Therefore extra-flat directions are essential for an MSSM flat direction to play 
the role of the inflaton in D-term inflation. 



E. Post-inflationary evolution 



Including soft SUSY breaking terms, the potential is 

Ai0 m 



Ai0 



H.c. 



Wt\4>- 



2d 



m—1 



M 



I m— 3 



2 2 



|2\2 



(26) 
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where we used equation (J6j). Here m| is the soft SUSY breaking scalar mass. The potential 
has two SUSY non-renormalisable terms: the third term g 2 £ 2 (|0_| 2 /£)|0 n /0™| 2 and the 
fourth term |Ai0 m /M m_3 | 2 . Assuming that m — 1 < n (since the case of most interest will 
be that where m = n), the third term is dominant if > 0*, where 

,n-m+i = MS ( £ \ ' (27) 

while the fourth term is dominant if < 0*. An important point in what follows is that 
for 0^0*, the v4-term will be effectively suppressed compared with the usual case of an 
MSSM flat direction with potential stabilised by a non-renormalisable term. This is because 
the A-term is coming from the first term in the superpotential, equation ([1]), whereas the 
non-renormalisable term in the scalar potential is from the second term. As a result, the 



baryon asymmetry generated by the Affleck-Dine mechanism [17J, [18] will be suppressed 
relative to the MSSM flat direction case. Once inflation ends, the 0_ field oscillates around 
the minimum (0_) = y/£ 6 . The field will oscillate around the origin dominated by either 
the |Ai0 m_1 /^ m ~ 3 | 2 or the | (0_) \ 2 g 2 ^\4> n / 4>^\ 2 term, depending on the amplitude. While 
the amplitude of the oscillation is large, the energy density of will decrease more rapidly 
than that of 0_ (V oc <p d implies that p oc a ^ 6d ^ d+2 \ with d > 6 for oscillations and d = 2 
for 0_ oscillations), so the Universe initially becomes 0_ dominated. In the following we will 
assume that the 0_ oscillations efficiently decay into radiation. (We will comment on how 
our results are altered if this is not satisfied.) Due to the 0_ decay, the radiation produces 
two distinct thermal corrections to the potential equation ( 1261) . The field is expected to 
acquire a thermal mass term 

/i 2 T 2 |0| 2 , (28) 



with h being a coupling between and a particle in the thermal bath [19(] in the case where 
the expectation value of the field is relatively small and the radiation temperature is high 
enough, and also a logarithmic term 

«T 4 ln^, (29) 
which appears at the two-loop level through the running of couplings with non-vanishing 



201 ] . Here, a is a constant of order of 10 2 and its sign can be positive or negative. For 



In general, the minimum of the potential is at = £ 1/2 (l - |0/0 c | 2n ) 1/2 . This rapidly tends to 

= £}/ 2 as the 4> oscillations are damped from 4>c to small amplitudes. 
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our purpose, hereafter we consider the case that a is negative. (For a positive a, the field 
oscillates around the origin by either the mass m^, the thermal mass or this two-loop effect 
and simply decays into radiation.) The potential with the two-loop induced logarithmic 
potential is 



V{<j>) 



1 + K\n 



A 2 



+ Am 



3/2 



Ai0 r ' 



mM 



m— 3 



+ H.c. + aT 4 ln 



J^2 



<j) n 


to 








M m ~ 3 



(30) 



where we include the radiative correction to with the direction dependent coefficient K. 
or MSSM flat directions which do not include a large top quark component, K ~ — 10~ 2 



21 



221 ] . However, a large top squark component can drive K to positive values. 



A negative K is the source of the spatial instability which leads to Q-ball formation in 
the gravity mediated SUSY breaking model 2l|. For a negative a, as shown in Ref. 23], 
the thermal mass term cannot appear because of a relatively large expectation value of 
the field. Here, the <fi field is trapped with nonvanishing value by the thermal logarithmic 
term, equation (1291) . and the non-renormalizable term, until the temperature decreases to a 
certain value. As the temperature falls, the expectation value of <j) becomes small. When 
the (t> becomes as small as 

(-a)T 4 

(31) 



cf) starts to oscillate around the origin with the angular momentum in the space induced 
by A-term, which is equivalent to the charge density (baryonic and/or leptonic) carried by 
4> |17|]. Provided that the reheating by the 0_ decay is completed before cf> starts to oscillate, 
we find from equation (|3T|) that 

:5(l, -"> (32) 



PR 



n 2 g* 



Here, is the effective total degrees of freedom of the relativistic species in the radiation. 
Since the ratio in equation ( 132]) is of order of 10 -4 , the field oscillations (or the Q-ball 
density formed from the condensate if K < 0) soon dominates the Universe. 



III. REHEATING AND BARYOGENESIS 



Reheating in this model is from the decay of the extra-flat direction inflaton field. The 
reheating temperature will therefore depend on whether or not the flat direction condensate 
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xag ments into Q-balls, which in turn depends on the t squark content of the flat direction 
221. 



A. Q-ball formation and decay 



The field oscillates around the origin coherently to begin with, but there is a spatial 
instability of its fluctuations due to the negative K. After inhomogeneities in the field grow, 
the coherent fragments and, as a result, Q-balls are formed 2jJ, |24j. Here we briefly 
summarize properties of Q-balls in gravity mediated SUSY breaking models. The radius of 
a Q-ball, R, is estimated as R 2 ~ 2/(\K\m 2 b ) [21|. By numerical calculations, it was shown 
that almost all the produced charge is stored inside Q-balls, and that a good fit to the Q-ball 
charge is 



m 



with 



for 



> 



for e < e r 



and 



rrir, 



2q\A\ ( ) suiA x Mm 
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Ai * 

6 x icr 3 



rrir, 



(33) 



(34) 



(35) 



25] . (For e < e c the 



where 5 is the CP violating phase, e c ~ and j3 
condensate will fragment to pairs of oppositely charged Q-balls.) The last two factors in 
equation ( |35l) are, respectively, the suppression of the baryon asymmetry due to the effective 
suppression of the A-term relative to the non-renormalisable term once 0^0* (where Ai * is 
defined below), and the enhancement due to H os <ti at the onset of oscillations, which 
allows the B violating A-term to act over many oscillations before expansion diminishes 
the A-term 7 . 

The decay temperature of Q-ball is given by 



26] 



V2 no 



VlTeV/ 



20 \ 1/2 



Q 



GeV, 



(36) 



7 In the case where the (j>- field does not decay efficiently to radiation, the Universe after inflation will 
be dominated by <f>- oscillations and onset of oscillations will be typically determined by an order H 2 
correction to the <f> mass squared due to non-minimal Kahler interactions of the form |0_| 2 |</>| 2 . In this 
case ra^ w H os in Eq. 
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where 10 3 > / s > 1 is the enhancement factor in the decay if Q-balls can decay into final 
states consisting purely of scalar particles. Since Q-balls come to dominate the Universe in 
our scenario, the decay temperature gives the reheating temperature at the onset of radiation 
dominated Universe. The resultant emitted charge to entropy ratio is given by 



71,, 



(37) 



B. Affleck-Dine baryogenesis 



The baryon asymmetry is generated by the B and CP violating A-term when <fr starts to 
oscillate around the origin 171 ] . As usual in D-term inflation, there is no order H correction to 
the A-term before starts oscillating, since 0+ = throughout 27j. In addition, in the case 
of a non-singlet inflaton there can be no linear coupling of the inflaton I to superpotential 



monomials W in the 
A-term correction 
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Mahler potential of the form VW ', which would generate an order H 
]. Therefore the phase of the inflaton relative to the A-term at the 
onset of (f) oscillations, 9, is determined by its initial random value during inflation, in which 
case sin S ~ (sin2#)/2. This phase gives the CP violating phase required for Affleck-Dine 
baryogenesis, with tib oc 9 for 9 small compared with 1. In the estimation of the resultant 
baryon asymmetry produced by the Affleck-Dine mechanism, the important quantity is the 
amplitude of the AD field when it starts to oscillate, OS . For OS > 0*, the amplitude is 
given by 



ml ~ 



2^2 



/,n-l 

OS 



(38) 



where we assume |0_| 2 = £. On the other hand, for 

Ai0; 



the amplitude is given by 



m-2 

OS 



M 



m— 3 



(39) 



The former applies in the case of a small Ai ^ Ai *, with 

m-2 



(n-1) 




) M (m-3)(n-l)_ 

n 



n — m + 1 
2 \ 2 



(40) 



while the latter corresponds to a large Ai ^ Ai *. For the case Ai ^ Ai * we obtain 



|0o 



|2n 



l 

n-l 



(41) 
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Equation (133]) then gives 



\2n 



1 

7-1-1 



Q^ L q2 . 2 2fr-2) 1 ^ ( 4 2) 



The expansion rate at the onset of oscillations during radiation domination can be 
obtained from equations (J3~I|) and (j4"Tj) . 

# os f n 2 g*\ 1/2 |0 OS | ^ 



(44) 



m V 90a / M 
with 

i 



For n = 4 this gives, 



Therefore H os /m^ ~ 10~ 3 . This justifies neglect of corrections to the soft SUSY breaking 
terms at the onset of (ft oscillations. 
For n — 4, equation (l4"2l) becomes 

2i /0.1\ 2/3 /4 x 10 15 GeV\ 4/3 / \(ft c \ \ 8/3 /lTeV\ 4/3 , . 



Then from equations (1361) . (13~Tj) and (J43j) . and using m,p/H os « 10 3 , the decay temperature 
and baryon asymmetry are given by 

-V3 />, v mis n„i/-\ -2/3 / I 1 I \ -4/3 



and 



- 1/3 '4 x 10 15 GeV^ ~ 2/3 ' ^ ' N ~ 4/3 



y \0.3M 

(48) 

The observed baryon asymmetry is n q /s = (1.8 ± 0.1) x 10~ 10 . Hence e < 10~ 7 is necessary 
to account for the observed B asymmetry. The Q-ball decay temperature, which gives the 
reheating temperature, is in the range 3-100 GeV for 1 < / s < 10 3 . For Ai ^ Ai*, from 
equation ( 1331) we have e ~ (0.1 — l)(Ai/Ai ^)(m^/ H os )9. The random phase of (ft during 
inflation would be expected to be of order 1; therefore in order to generate the observed 
baryon asymmetry we require that Ax ~ (10~ 10 — 10~ 9 )Ai *. With m = n = 4, Ai * is, 

, lt ,,3x 10 - 8 (^)-(i2~)- 2/ \f. (48, 
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Since A2 should not be very small if g is not very small, Ai * ~ 10~ 8 is likely. Therefore 
to account for the observed baryon asymmetry with 9 pa 1 we must have Ai < 10 -17 . This 
corresponds to suppression of the QQQL or u c u c d c e c superpotential terms by a mass scale 
A > 10 35 GeV. Thus a much larger suppression is necessary for successful baryogenesis 
than is required by proton stability Note that it may be possible for Q-balls to decay at 
a temperature greater than that of the electroweak transition if / s pa 10 3 , corresponding to 
Q-ball decay to purely scalar final states. In this case any dangerous baryon asymmetry will 
be erased by B + L violating sphaleron fluctuations. 



C. Baryon isocurvature perturbations 

The CP violating phase 5 is given by the phase of <ft during inflation relative to the 
A-term, which defines the real direction. Therefore in the case where reheating is via Q- 
ball decay, e ~ 10 2 (Ai/Ai *)6 implies that 9 pa 1CT 9 (Ai */Ai) in order to have e ~ 1CT 7 , as 
required for successful baryogenesis with / s = 1. Quantum fluctuations of in the phase 
direction will lead to baryon isocurvature perturbations, which can be large when < 1. 
For uncorrelated baryon isocurvature perturbations, the fractional contribution to the CMB 



power spectrum is given by «bi, where 28|, |29j 



" BI = fej ■ (50) 

with 5n B /riB ~ fo59. In our case fg ~ 1/9. The present observational limit is «bi < 0.26 
fl- With 6 w 10- 9 (Ai ,/Ai), fi DM = 0.23 and Vl B = 0.04, this gives a upper bound on H/<f>, 

H < 1n _ia f^i 



^ 10- 13 ^ . (51) 

t^ 1 / 2 ~ 4 x 10 15 GeV, corresponding to O(10)% cosmic strings, implies that H = 2.7 x 
10 12 ^ GeV. Since we are considering pa <p c m (0.1 — 0.3)M, we therefore have 

pa (0.6 - 1.8) x 1(T 6 2 , (52) 

Thus with g pa 0.01 — 0.03, the baryon isocurvature perturbation is sufficiently small if 
A1/A1 * ^ 10~ 5 . The correct baryon asymmetry then requires that 9 ^ 10~ 4 . Thus even if the 
initial random phase of the flat direction field could satisfy 9 <C 1, the flat direction would 
still have to be suppressed by A > 10 31 GeV in order to avoid large baryon isocurvature 
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perturbations. For A pa 10 31 GeV and 9 pa 10 4 , the correct baryon asymmetry will be 
generated together with a potentially observable baryon isocurvature perturbation. 

D. Non-thermal dark matter 

The reheating temperature is ~ 1 GeV for the case where Q-ball decay to purely scalar 
final states is kinematically suppressed, such that / s = 1. This low reheating temperature 
implies that Q-balls may decay below the freeze-out temperature of neutralino LSPs, in 
which case Q-ball decay will also produce non-thermal LSP dark matter particles. In fact 
dark matter particles are often overproduced, in particular for the standard bino-like neu- 
tralino LSP. Although several way s to avoid this problem have been proposed by taking an 
alternative choice of the LSP 23j, |3l|, |32|, |33j, perhaps the simplest ones are to assume a 
Higgsino-like neutralino LSP [sj, or a gravitino LSP [3] with a sneutrino NLSP to escape 



BBN constraints 34]. 



E. Reheating from flat direction condensate decay without Q-ha\l formation 

In the case where the inflaton corresponds to a flat direction with a large t squark compo- 
nent, the <p condensate will not fragment to Q-balls since K > [22j. In this case reheating 
will occur via conventional flat direction condensate decay and a higher reheating temper- 
ature is expected. For the B — L conserving u c u c d c e c and QQQL directions, the baryon 
asymmetry from Affleck-Dine baryogenesis will be erased by sphaleron B + L violation so 
long as the condensate decays at T > T ew . In this case it is possible for the initial phase 
of (j) to take its natural value, 6 ~ 1, without requiring a suppression of the flat direction 
beyond that required to evade proton decay. Assuming that <fi oscillations dominate the en- 
ergy density when the <fi field decays to radiation, the energy density is given by p pa m\<f>\, 
where is the amplitude of the oscillations when they decay. \<fi&\ is then related to the 
decay temperature T d by 



2 _ fc d T d 4 _ 7r 2 g(T d 

ml 30 



= • (53) 



For h\(j)d\ > rn^, where h is the gauge or Yukawa coupling of MSSM particles to the flat 
direction, particles coupling to <fi gain masses greater than and so the <fi decay is kinemat- 
ically suppressed. Therefore the condensate decays once ~ m^/h, assuming that T d > H 
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when this occurs. The energy density in the field at this time is p ps m^/h 2 . Therefore the 
decay temperature, which is equivalent to the reheating temperature Tr, is 

T '~~^r- (54) 

With g(T D ) ps 200 we find k d ps 65. Therefore 

where the particles with the smallest coupling h to <p will dominate the decay process, so 
long as Td > H . Therefore T d ~ 1 — 100 TeV in this model, assuming that the smallest 
coupling satisfies 0.1 > h ^ 10~ 5 . Once h<f> < the decay rate may be estimated to be 
Ta ~ /i 2 m^/47r, so the condition > H ps 5T 2 /M is easily satisfied for in this range. 

We have assumed that the kinematic suppression of the decay rate prevents <fi decaying 
until > h(p, in which case T<j ^ 100 TeV. We should check that decay through heavy 
intermediate particles cannot cause it to decay significantly earlier. The decay rate via heavy 
intermediate particles of mass tup will have the generic form 

Td w W ' (56) 

where «d < 1 is a product of couplings and phase space factors. Since there are two heavy 
intermediate states, r > 4 is expected. For r = 4, and using equation (}53"j) . (J56"j) and 
T d ~ H(T d ) ps BT^/M, this gives for the decay temperature 

/ \ 1/10 / \ V 10 , s Q/m 

^ - ( -^"1 MM) 1 ' 10 - 35 f -^O TeV . (57) 

\g A k 2 d k T J y * > Wk 2 d k T J vl TeV/ 1 ; 

Thus for typical couplings, the decay through intermediate states will also result in a re- 
heating temperature in the range Tr ~ 1 — 100 TeV. It is significant that the reheat- 
ing temperature, Tr < 100 TeV, is naturally compatible with the thermal gravitino upper 
bound, Tr < 10 6 GeV, without any tuning of couplings. Even though the inflaton is part of 
the MSSM sector, it still leads to the required low reheating temperature. Since sphaleron 
B + L violation will erase the baryon asymmetry produced by the flat direction inflaton de- 
cay, baryogenesis must occur via some other mechanism, such as Affleck-Dine baryogenesis 
along an orthogonal flat direction. 
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IV. CONCLUSIONS 



We have shown that it is possible for an MSSM extra-flat direction (one suppressed by an 
effective mass scale much larger than the Planck mass) of the form QQQL or u c u c d c e c to play 
the role of the inflaton in a D-term inflation model. This eliminates the otherwise unmo- 
tivated singlet inflaton, reducing the number of required additional fields and so providing 
a minimal version of D-term inflation. The model has all the advantages of conventional 
D-term inflation with respect to compatibility with SUGRA and absence of fine-tuned cou- 
plings. The nature of reheating depends on whether the extra-flat direction is unstable with 
respect to Q-ball formation. In the case where Q-balls form, it is possible to generate the 
baryon asymmetry via Q-ball decay so long as the mass scale suppressing the flat direction 
is sufficiently large, A > 10 31 — 10 35 GeV, depending on the random phase 6 of the flat direc- 
tion scalar during inflation. With A « 10 31 GeV and 9 m 10~ 4 it is possible to generate an 
observably large baryon isocurvature perturbation. The reheating temperature from Q-ball 
decay is typically in the range 3 — 100 GeV. As this can be less than the neutralino LSP 
freeze-out temperature, it is also possible to produce non-thermal dark matter from Q-ball 
decay. In the case where the flat direction has a large t squark component, there is no Q- 
ball formation. In this case the reheating temperature from decay of the homogeneous flat 
direction condensate is in the range 1 — 100 TeV, ensuring sphaleron erasure of the baryon 
asymmetry from the B — L conserving directions while remaining naturally compatible with 
the thermal gravitino upper bound on Tr. The fact that we are able to calculate the re- 
heating temperature in this case is a direct consequence of the inflation being part of the 
MSSM sector. Since the baryon asymmetry from the flat direction is erased, the mass scale 
suppressing the flat direction in this case is constrained only by proton decay, A > 10 26 GeV. 
We have interpreted the WMAP observation of the spectral index as being due to an order 
10% CMB contribution from cosmic strings combined with a nearly scale-invariant adiabatic 
curvature perturbation, n s ~ 1. As in conventional D-term inflation, we can simultaneously 
suppress the contribution of the cosmic strings to the required level while increasing n s by 
considering a small enough coupling of the inflaton to the Fayet-Iliopoulos charged fields. 
The extra-flat direction D-term inflation model has an advantage over conventional D-term 
inflation in that the range of coupling which leads to an order 10% contribution from cosmic 
strings is much wider, making it perhaps more natural. For this solution to work, it is also 
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necessary to have a U(l) gauge coupling that is somewhat smaller than the known gauge 
couplings, g ~ 0.01 — 0.03. A significant feature of the model is that the superpotential 
coupling of the monomial QQQL or u c u c d c e c to $ + $_ must be much larger than the pure 
monomial superpotential coupling. This feature may serve to test the compatibility of the 
model with an ultra-violet complete theory, as we would naively expect all the superpotential 
couplings of the monomial to be strongly suppressed. Finally, we note that other solutions to 
the cosmic string and spectral index problems are possible, for example SUGRA corrections 
from a non-minimal Kahler potential 35] and/or modification of the inflaton potential by 
other fields, such as a RH sneutrino 36| . 

The existence of extra-flat directions of the MSSM is one way to interpret the empirical 
suppression of non-renormalisable MSSM superpotential terms demanded by proton stabil- 
ity. It will be important to establish whether extra-flat directions can be understood in 
the context of an ultra-violet complete theory and to explore more generally their role and 
possible signatures in cosmology. 
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